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We critically examine the Hartree-Fock-Bogoliubov (HFB) solution of the equations of motion 
for condensate fluctuations in a weakly coupled plasma of charged bosons at finite temperature. 
Analytic expressions are derived for the first two infrared-divergent terms in both the momentum 
distribution of the noncondensate and the anomalous Bose correlation function at low momenta, 
fncorporation into the theory of the appropriate form of the Hugenholtz-Pines relation for the 
chemical potential is needed to cancel an unphysical divergence. Exact cancellation of infrared- 
divergent terms is demonstrated in the HFB shift of the single-particle excitation energy away from 
the Bogoliubov value at long wavelengths, with the residual terms raising it towards the plasma 
frequency at low temperature. Numerical illustrations are presented for a number of properties of 
the boson plasma as functions of temperature and density in the weak-coupling regime: these are 
the chemical potential, the condensate fraction, the normal and anomalous momentum distribution 
functions and the corresponding one-body density matrices, and the dispersion relation of single- 
particle excitations. 

PACS numbers: 05.30.Jp, 03.75.Fi, 74.20.-z 



I. INTRODUCTION 

The fluid of charged bosons over a uniform neutralizing background (CBF) interacting via inverse-first-power 
repulsions is a model in quantum statistical mechanics! having possible relevance to superconductivity in layered 
cuprate materials! and to the equation of state and nuclear reactions in dense plasmas of pressure-ionized helium.! Its 
counterpart with inverse-first-power attractions can arise in systems of atoms immersed in particular configurations 
of intense off-resonant laser beams.! 

The properties of the CBF at zero temperature have been addressed in a number of theoretical studies, ranging from 
the early work on the weak-coupling regimcl to treatments based either on the use of Jastrow-Feenberg correlated 
wave functions! or on self-consistent accounts of correlations.! Quantum Monte Carlo calculations! have covered the 
whole range of coupling strength up to Wigner crystallization. More recently, the evaluation of the condensate fraction 
has been extended to finite temperature! within the theoretical frame provided by the Bogoliubov formalism. These 
theoretical results could only be tested against Monte Carlo data at zero temperature, showing that this most simple 
approach yields almost quantitative values for the interaction-induced depletion of the condensate over a non- negligible 
range of coupling strength. 

However, the standard Bogoliubov theory can be expected to fail with increasing temperature, as the noncondensate 
fraction increases. The purpose of the present work is to examine the roles of both the noncondensate and the anoma- 
lous Bose correlations in the CBF at finite temperature within the Hartee-Fock-Bogoliubov (HFB) approximation, 
with due attention to the chemical potential of the fluid as given by the Hugenholtz-Pines relation.EJ 

It may be mentioned at this point that interest in the HFB approach to boson fluids at finite temperature has been 
revived in relation to developments in the study of Bose-Einstein condensed atomic gases under confinement. In this 
context GriffirB has critically analyzed the HFB scheme with special regard to the classification of (conserving versus 
gapless) approximations given by Hohenberg and MartinEl for fluids of neutral bosons. The HFB approach has also 
been used to describe the shape deformation modes of confined atomic gases in dependence of temperature. 

In his work Griffin0 has emphasized that the problem of determining the excitation spectrum of a Bose-condensed 
fluid is difficult because of the need to satisfy several different constraints. In particular, in an explicit treatment of a 
Bose gas of neutral particles with contact interactions he has stressed that the HFB value for the chemical potential 
violates the exact Hugenholtz-Pines expression for this quantity and that the predicted excitation spectrum violates 
the gaplessness requirement. Proportionality between the single-particle and the collective sound-wave excitation 
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spectra in the neutral fluid at long wavelengths can be understood from the argument developed by Gavoret and 
Nozieres0 The long-range nature of the Coulomb interactions in the CBF affects these properties and emphasizes 
the shortcomings of the HFB approach. Because of the plasmon gap in the collective excitation spectrum at long 
wavelengths, one expects the same gap in the single-particle spectrum, at least at T = 0.00 As we shall sec, the 
violation of the Hugenholtz-Pines relation is associated in the CBF to infrared divergencies and requires ad hoc mend- 
ing. Analytical results can then be derived in the long-wavelength (low-momentum) region, extending exact results 
previously obtained^ from sum-rule arguments at zero temperature. The modified HFB approach, incorporating the 
Hugenholtz-Pines value of the chemical potential, is amenable to numerical evaluation and in particular displays the 
expected gap in the single-particle spectrum of the CBF at low temperature, though a discrepancy from the value of 
the plasmon gap remains. 

In brief, the plan of the paper is as follows. Section II collects for convenience the main standard equations of 
the HFB approach and presents their modification that we propose for a Bose plasma in order to incorporate the 
Hugenholtz-Pines relation. A further subsection presents our analytical results for the long-wavelength region and 
our numerical results for the limiting value of the excitation energy as a function of temperature in the same region. 
Various other numerical results are presented in Sec. Ill, while Sec. IV gives a brief summary and our conclusions. 



II. THE HFB APPROACH TO A BOSON PLASMA 

As is well known,0 the ideal Bose gas starts undergoing Bose-Einstein condensation at a critical temperature 
), where n is the particle number density and m the particle mass. The condensate density no 
increases with decreasing temperature, to become equal to the total density n in the ideal gas at zero temperature. 
The role of the interactions is to depress the condensate fraction and to shift the critical temperature T c away from 
To. Within the HFB approach this shift is a coupling-strength-dependent increase in the present case of Coulomb 
repulsions. 

The CBF is described by the Hamiltonian 

H = J dr^{r) (-|^V 2 - f?j ip{r) + \J J drdr'^(r')^(r)V(\r - r'\)ip(r')ip(T) (1) 

where ip(r) is the bosonic field operator, /i is the chemical potential and V(r) is the Coulomb interaction potential. 
The average potential felt by each particle vanishes because of the presence of the neutralizing background: that is, 
in the calculations reported below we take the Fourier transform of V(r) as V4 = Aire 2 /k 2 with Vk=o set equal to 
zer 00. The coupling strength is measured by the dimensionless parameter r s , defined by r s as — (4am/3) with 
as the Bohr radius. 

The HFB approach is based on the equation of motion for the field operator and adopts a decoupling procedure in 
order to linearize it and solve it. The Bogoliubov prescription starts by separating out the condensate part from the 
field operator ip(r) by setting V>( r ) = "4>o a o + </>( r )7 where So and 4>(r) are the field operators for the condensate and 
for the cloud of particles promoted out of the condensate. Here, ipo = (VK 1- )) is the condensate wave function, given 
by the ensemble average of ijj(r) in the symmetry-broken state. In the homogeneous gas we take the condensate wave 
function as uniform in space and time, and for the sake of simplicity we also assume that it is real. 

The equation of motion obeyed by the Heisenberg field operator is0 

f_|lv 2 - u W)+ / dr'^(r',t)i;(r',t)^(r,t)V(\r-r'\). (2) 



dt \ 2m 

Following a standard procedure, we use the Bogoliubov prescription in Eq. (2) to derive the equation of motion for 
the fluctuation operator <fi(r). 

In the HFB scheme the following " mean field" approximation is adopted,Q 
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<£t( r ',^(r,i) S (0+(r',t)0(r,t)) = rz(|r - r'|,t) 
^(r', t)0(r, t) <* (0(r', t)0(r, i)) = m(|r - r'|, t) 



(3) 



where we have introduced the normal and the anomalous density matrix, h and m. The three-point kernel 
ffl(r', t)(j>(r', t)cf>(r, t) is approximated as 

0V, *)<H r ', *)<X r , *) = *)<K r ', *)) <^ r , *) + ri(|r - r'|, t)4>(r', t) + m(|r - r'|, t)^( r ', t ) (4) 

The equation of motion for the fluctuation operator in the HFB approximation then reads 

th ^A = - M ) 0(r,t) + | rfr' [n + n(|r - r'|,i)] 0(r',t) + [n + m(|r - r'|,t)] ^(r'^Mk - r'|). (5) 

By taking the average of Eq. (2) and using the approximations set out in Eqs. (3) and (4) we also obtain the 
expression for the chemical potential in the HFB approximation, 

[Ihfb= J dv'[n(\r-r'\,t) + m(\r-v'\,t)]V(\r-r'\). (6) 

The condition of charge neutrality has been used in deriving Eqs. (5) and (6) by setting V k= a — 0, as already 
remarked. 

The solution of Eq. (5) is achieved by means of a normal-mode expansion for the fluctuation operator, 

4>(r, t) = [uj(r) cxp (-iEjt/h) aj + v*(r) exp (iEjt/K) a] . (7) 
3 

This leads to a set of Bogoliubov - de Gennes equations, which can be diagonalized in Fourier transform to obtain 
the expression for the HFB spectrum of single-particle excitations, 



E k = yJ[e k + I n (k)] 2 -I m (k) 2 (8) 



and for the HFB mode amplitudes, 



u 2 = \[l + {e k + I n {k))/E k ] 
v 2 = \[-l + {e k + I n {k))/E k ] 



(9) 



In these equations e k = (h 2 k 2 /2m) — fi and the self-energies {/„(&), I m (k)} are to be determined self-consistently 
with Eqs. (8) and (9) according to 

{I n (k),I m (k)} = n V k + Vj q _ k , {n(q), m(q)} (10) 

Here. 

n{q) = v 2 q + f q (u 2 + v 2 ) 
m(q) = u q v q (2f q + 1) 



(11) 



with fq = [exp (E q /kBT) — 1] 1 are the Fourier transform of the normal and the anomalous density matrix. Evidently, 
n(q) is the momentum distribution of the Bose fluid. Thus, the condensate fraction is given by 

n =n-J2[v 2 q +f q (u 2 + v 2 )] (12) 

q#0 
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(13) 



Finally, the normal and anomalous one-body density matrices are given by 

n(r) = n + £ q#0 [ v g + U { u \ + v g)] ex P ( ic l ' r ) 
m(r) =n Q + J2q^o K w <? ( 2 fq + *)] ex P («1 ' r ) 

Equations (8) - (10) reduce to the standard Bogoliubov approximation (see, e.g. Ref. [9]) when the contributions 
coming from the density matrices n(q) and m(q) are omitted. In the same approximation the chemical potential can 
be set to zero. In the HFB approximation, however, Eq. (6) for the chemical potential of the CBF is plagued by an 
infrared divergence (see below). 

We overcome this difficulty in the present case of a Bose plasma by recourse to the exact expression given for the 
chemical potential by the Hugenholtz-Pines relation. This states that 



H = En(q = 0, u = 0) - Eia(q = 0, u = 0) 



(14) 



where £jj(q, u>) is the 2x2 matrix for the self-energies of a Bose-condensed gas.E 2 ] We then replace the divergent 
expression in Eq. (6) by the expression 



(15) 



Hhp= J dr' [n(|r-r'|,i)-m(|r-r'|,i)]F(|r-r'|) 

q#0 



Equation (15) follows from Eq. (14) by using the HFB expression for the self-energies given in Eq. (10), together 
with the requirement of charge neutrality which suppresses the q = term in the sum. 

Equations (8) - (11) are well behaved when one adopts the expression in Eq. (15) for the chemical potential. Their 
solution will be presented in Sec. Ill below. However, before proceeding to numerical calculations we need to evaluate 
analytically the long-wavelength behavior of the CBF in the HFB approach. 

A. Long-wavelength behaviors 



The following asymptotic expressions will be demonstrated below to hold at low momenta for the CBF at temper- 
ature T > in the HFB approximation: 



lim \I n (k) 

k^O 



I m {k)} 



4ire 2 n (47re 2 



(2/o + 1) 



(16) 



fc 2 16E k 

where we have introduced the notations Eq — Ek=o and /q = fk=o- The infrared divergence in the HFB expression 
for the chemical potential in Eq. (6), and the cancellation of the divergence at both leading and subleading order in 
Eq. (15) from the Hugenholtz-Pines relation, are immediately evident. 

The corresponding asymptotic expressions for the normal and anomalous distribution are obtained from Eqs. (9) - 
(11), which yield 



lim \n(k) 
k^o 



m(k)} = (2/ 



1) lim [I n (k)/2E k ] 

k— >0 



The result is 



lim \n(k) 

k^O 



m{k)] = (2/ + 1) 



47re 2 n (47re 2 



n 



2£ fc 2 



32E%k 



(2/ + 1; 



(17) 



(18) 



This asymptotic expression for n(k) reduces at T = to the exact result obtained by sum rule arguments in Ref. [15], 
provided that the energy Eq coincides with the plasma frequency io v = W 1 Amne 1 jm. We present at the end of this 
section a calculation of Eq in the HFB approach and compare the result with the plasma frequency. 
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To demonstrate the results given above in Eqs. (16) and (18), we examine the HFB expressions for the self-energies 
{I n (k), I m (k)} in Eq. (10). The integral on the RHS of this equation is easily shown not to contribute at order k~ 2 , 
so that the leading-order term in Eq. (16) is immediately established. The leading-order term in Eq. (18) follows 
fromEq. (17). 

With the leading-order expression for n(k) and m{k) we can now return to the evaluation of the integral on the 
RHS of Eq. (10) and determine the subleading term in Eq. (16). For the quantity I n (k) we write 



X) Vj k _ q |n(g) =4ne 2 J2 



1 



q#0 



q#0 



k 2 - 2k ■ q + q 2 



n(q) — lim n(q) 



( W) 2 n n V 1 

2E Q [ J ° ' ^ Q q 2 (k 2 - 2k -q + g 2 ) 



(19) 



The value of the second integral on the RHS of Eq. (19) is l/(8fc), while the first integral does not contribute to this 
order.El With an identical argument for I m (k) the subleading-order terms in Eq. (16), and hence in Eq. (18), are 
thus established. 

We turn next to the evaluation of the energy Eq from Eq. (8) in the long- wavelength limit. We get 



Eo/h 



( 4ne 2 n 



1/2 



2m 

Inn 1 + -g — 



[I n (k) - I m (k) - m] 



1/2 



\ m 

We also have from Eqs. (10) and (14) that 

lim [I n (k) - I m (k) -fx} = \k 2 Hi) ~ m (l)} 



(20) 



(21) 



q#0 



where according to Eq. (11) we have 



»fe)-"(?) = ~ + |£(2/. + l) 



so that 



lim [n(q) — m(q)] 

q^O 



(22) 



(23) 



We can thus write the final expression for Eq in a form which is suitable for numerical calculation, 

1/2 



Eo/h 



1/2 



2m \ - n 



3h" 



q^O 



n(q) — m(q) + 



(24) 



This expression shows that the value of E is finite, as implicitly assumed in the foregoing derivation of Eqs. (16) 
and (18). The standard Bogoliubov approximation yields Eo/h = (47re 2 no/m) 1 / 2 . As is shown in the left panel in 
Figure 1, the Hartree-Fock terms inside the brackets in Eq. (24) raise the value of Eo/h towards the plasma frequency 
(47re 2 n/m) 1 / 2 at low temperature, without quite reaching it. In both the Bogoliubov and the HFB approximation 
the value of Eo/fi drops with increasing temperature to vanish at the critical temperature T c . The corresponding 
behaviour of the chemical potential from Eq. ( |l5| ) is shown in the right panel in Figure 1. 

It is also worth making a remark on the asymptotic result given in Eq. (23). The quantity n{q) — m{q) is the 
HFB value for a special component of the two-body density matrix in the CBF, that is n(q,q)/no- It was shown in 
Ref. [15] from an exact sum-rule argument that this quantity takes the value —1/2 at low momenta in the CBF at 
T = 0. Thus, this exact asymptotic result holds in the HFB approximation, which yields the same limiting value at 
all temperatures T < T c . 
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III. NUMERICAL RESULTS 



We present in this section the results of numerical calculations of some other properties of the weakly coupled CBF 
in the modified HFB approximation. Figure 2 reports the condensate fraction no/n as a function of temperature at 
coupling strength r s = 1. In comparison with the result of the simple Bogoliubov (B) approach, there is a further 
large shift in the value of the critical temperature for the formation of a condensate away from the ideal-gas value. 
At T = the change in the interaction-induced depletion of the condensate from the result of the simple Bogoliubov 
(B) approach is quite small at low r s (r s < 5, say), thus confirming the results reported in Ref. [9] and ensuring 
good agreement with the Monte Carlo data obtained by Moroni et aln However, with further increase of r s the HFB 
condensate fraction drops very rapidly, in disagreement with the Monte Carlo data. 

In view of the asymptotic behaviors given in Eq. (18), Figure 3 reports the quantity q 2 n(q) at T = in order to 
display the subleading term which is absent in the simple Bogoliubov approximation. This term is also evident in the 
Monte Carlo data, which are reported in Figure 3 from the work of Moroni et alB, and there is in fact reasonable 
quantitative agreement between the HFB momentum distribution and these data. The HFB predictions for the normal 
one-body density matrix n(r) at various temperatures are shown in Figure 4. 

Similar results are reported for the anomalous Bose correlations in Figures 5 and 6. Figure 5 shows both q 2 n(q) and 
q 2 m{q) at r s = 1 and two temperatures, to illustrate how departures arise with increasing momentum away from their 
common limiting values in Eq. (18). The normal and anomalous one-body density matrices are given in Figure 6. 

Finally, Figure 7 reports the HFB dispersion relation of single-particle excitations as a function of temperature at 
r s = 1. The deviations from the plasma frequency in the long- wavelength limit and their dependence on temperature 
have already been displayed in Figure 1. The Figure shows that, aside for this, there is very little dependence of the 
dispersion relation on temperature. There is, of course, no damping of these excitations within the HFB approach. 

IV. CONCLUSIONS 

Summarizing, we have examined the results given by the Hartree-Fock-Bogoliubov approximation for the quasi- 
particle excitations in a weakly coupled boson plasma through the temperature range in which a Bose condensate is 
present. The long-range nature of the Coulomb interactions emphasizes the difficulties of the theory that in previous 
work on neutral boson fluids have been shown to arise from subtle dynamical correlations induced by the presence of 
a condensate. For the boson plasma the Coulomb interactions generate infrared divergencies in both the normal and 
the anomalous momentum distribution, which are reproduced in the HFB approach at both leading and subleading 
order. However, the Coulomb case requires an ad hoc adjustment of the chemical potential. Exact cancellation of 
divergent terms occurs between the two distributions in the evaluation of the quasiparticle excitation energy at long 
wavelengths. The residual non-divergent HFB terms raise the quasiparticle energy towards the plasma frequency at 
low temperature, without quite satisfying the basic requirement that the Coulomb-induced gap in the quasiparticle 
spectrum at long wavelengths should coincide with the plasmon gap in the collective density-fluctuation spectrum. 

It appears, therefore, that the HFB approach may not provide an entirely fruitful starting point for an internally 
consistent description of the boson plasma. An approach based on the so-called dielectric formalismEl may be more 
promising and we hope to examine it in future work, with specific attention to the relationship between collective and 
single-particle excitation frequencies as a function of temperature. 
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T/T c T/T c 
FIG. 1. Left: the HFB single-particle excitation frequency Eq at long wavelengths (in units of the plasma frequency ui p ) as 
a function of temperature T/T c (in all figures T c is the HFB critical temperature) for three values of the coupling strength r s . 
The dots marked by B show the corresponding values in the standard Bogoliubov approximation at T = 0. Right: cutves for 
the chemical potential from the Hugenholtz-Pines relation. The QMC values reported at T = are from Moroni et alu 
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FIG. 2. The HFB condensate fraction no/n at r B = 1 as a function of temperature T (in units of e 2 /asA;B)), compared with 
the results for the ideal boson gas (Free) and for the Bogoliubov approach (B). 



8 




9 



10 



qa B 



FIG. 5. The HFB anomalous momentum distribution q 2 m(q) at r a — 1 and two values of T/T c , compared with the HFB 
normal momentum distribution q 2 n(q). Both curves are in units of a fl 2 . 
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FIG. 6. The HFB normal and anomalous density matrices at r s — 1 and two values of T/T c . 
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FIG. 7. The HFB dispersion relation for single-particle excitations (in units of the plasma frequency co p ) at r s = 1 and 
various values of T/T c . 
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